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Abstract 

Two-pion Hanbury-Brown-Twiss (HBT) correlations for central Pb-Pb collisions at the Large- 
Hadron-Collider (LHC) energy of y/s^m = 2.76 TeV are investigated for the first time with the mi- 
croscopic transport model UrQMD (Ultra-relativistic Quantum Molecular Dynamics). The trans- 
verse momentum dependence of the Pratt-Bertsch HBT radii is extracted from a three dimensional 
Gaussian fit to the correlator in the longitudinal co-moving system (LCMS). Qualitative agreement 
with the ALICE data is obtained, however Rout is overpredicted by nearly 50%. The LHC results 
are also compared to data from the STAR experiment at RHIC. For both energies we find that 
the calculated Rq/Rs ratio is always larger than data, indicating that the emission in the model 
is less explosive than observed in the data. 

PACS numbers: 25.75.-q, 24.10.Lx, 25.75.Gz 



I. INTRODUCTION 

In order to shed light on a large number of unsolved questions in fundamental physics 
l|-l3|, the Large Hadron Collider (LHC) at CERN had been designed, installed, tested, and 
repaired in the past two decades and finally, started normal operation in the end of the 
year 2009. Since then, a tremendous amount of experimental data in various aspects of 
high energy physics has been obtained and received much attention by theoretical physi- 
cists. Although the most exciting prediction, the Higgs boson, has not been observed in the 
data from the LHC experiments ^, the extracted bulk properties of the high temperature 
fireball created in such ultra-relativistic collisions have provided unprecedented information 
for fundamental investigations of the phase diagram of Quantum Chromodynamics (QCD). 
Here we want to explore the expansion properties of the created matter by investigating the 
spatial shape of the fireball. Although it is known that one can not measure the emission 
time pattern and/or the spatial profile of the source directly, a well-established technique, 
called "femtoscopy" or "HBT" (see e.g. 5| and references therein) can be employed to ob- 
tain this information. Femtoscopy has been extensively used in the heavy ion community 
since it provides the most direct link to the lifetime and size of nuclear reactions. The AL- 
ICE collaboration has published first results of two-pion Bose-Einstein correlations in both 
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p-p pI and central Pb-Pb [7J collisions at LHC energies in the beginning of the year 2011. 
These experimental results have attracted the research interest of several theoretical groups 
[THlOj, whose models are based on hydrodynamic/hydrokinetic and hadronic microscopic 
approaches. 

In this paper we show for the first time results for the HBT radii of two-pion correlations 
from central (< 5% of the total cross section cxt) Pb-Pb collisions at the LHC energy 
y^J/viv = 2.76 TeV from the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) 



model llNl4|. The calculations are compared to ALICE data as well as to those at the 



RHIC energy ^/sNN = 200 GeV. The UrQMD calculation results for p-p collisions at LHC 



energies are presented in |l5|. 



The paper is arranged as follows: In Section 2, a brief description of the UrQMD model 
and the treatment of the HBT correlations as well as the corresponding Gaussian fitting 
procedure is shown. Section 3 gives the main results of the model calculations. Finally, in 
Section 4, a summary is given. 



II. BRIEF DESCRIPTION OF THE URQMD MODEL AND THE HBT GAUS- 
SIAN FITTING PROCEDURE 

UrQMD |13l . Il4| is a microscopic many-body approach to p-p, p-A, and A-A interactions 
at energies ranging from SIS up to LHC It is based on the covariant propagation of mesons 
and baryons. Furthermore it includes rescattering of particles, the excitation and fragmen- 
tation of color strings, and the formation and decay of hadronic resonances. At LHC, the 
inclusion of hard partonic interactions in the initial stage is important and is treated via the 
PYTHIA pj model. 

In the present study, the cascade mode of the latest version (v3.3) of UrQMD is used (for 

nn 

details of version 3.3. see [13|, llJ]). Some predictions and comparison works with data from 
reactions at LHC have already been pursued based on this version and showed encouraging 
results for the bulk properties [l7|, Il8 |. 

To obtain HBT radii, first, about 200 and 10000 central events are calculated for Pb-Pb 
collisions at LHC and for Au-Au at RHIC, respectively. Then, all particles with their phase 
space distributions at their respective freeze-out time (last collisions) from UrQMD are put 
into an analyzing program using the formalism of the well known "correlation after-burner" 



(CRAB) 19|. At last, the constructed two-pion HBT correlator (regardless of charge) in the 
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without influence of residual interactions 



longitudinally co-moving system (LCMS) 

is fitted (using the x^ method) with a three-dimensional Gaussian form expressed as 



C(q,K) = l + A(K)exp 



(1) 



i,j=0,S,L 

In Eq. ([1]), A represents the fraction of correlated pairs [5l| and g^ is the pair relative mo- 
mentum q = pi — p2 in the i direction, pi being the momenta of the particles. L is the 
longitudinal direction along the beam axis, O the outward direction along the transverse 
component of the average momentum K of two particles (k^ = |pit + P2t|/2) and S the 
sideward direction perpendicular to the afore mentioned directions. The effect of cross terms 
with i ^ j on the HBT radii is found to be negligible in the present fits when a pseudorapidity 
cut \ri\ < 0.8 is used, as in experiments, and is not discussed in this paper. 

For central collisions, the HBT radii are, except for an implicit Kt dependence, related 
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FIG. 1: Projections of the three-dimensional correlation function (points) and of the respective 
fit (lines) for the fcr bin 200 — 300 MeV/c and \r]\ < 0.8. When projecting on one axis the other 
two components are restricted to the range (-30 < q < 30) MeV/c. 
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where x, y, z and t are the spatio-temporal separation of the particles in a pair and (3= 
K/i^o- If no space-momentum correlations are present the regions of homogeneity and the 
source size coincide. In central collisions the relation (x^) ~ (y^) is satisfied. Thus -Rq and 
i?| differ mainly in the last two terms of Eq. (jj]) . The first of these two terms is dependent 
on the strength of the correlation of emission time and transverse emission position, while 
the second one is especially sensitive to the particle emission duration. 

III. RESULTS 



The correlation functions are studied in bins of the transverse momentum kx = |kr|- 
Fig. [1] shows the projections of the three-dimensional correlation function (points) and of 
the respective fit (lines) for the kx bin 200 — 300 MeV/c. It is seen clearly that the correlator 
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FIG. 2: (Color online) (a) the transverse position Xout versus emission time distribution of pions 
with cuts \r]\ < 0.8, 200 < kr < 300 MeV/c, \qi\ < 100 MeV/c, and tfreezeout < 199 fm/c. (b) the 
contribution of terms < x'^ > (line with crosses), —2 < fSxtx > (line with asterisks), and < /3^t^ > 
(line with open circles) in Eq. (j2]), to Rq (line with triangles) with same cuts as for (a). 

in sideward direction can be described by a Gaussian form fairly well. However, it deviates 
slightly from a Gaussian in the other two directions, especially in the longitudinal direction, 
as found and discussed in previous publications for HICs at lower energies 23|. At LHC, the 
fraction of excited unstable particles is much larger than at lower beam energies, therefore, 
the non-Gaussian effect is more pronounced in the current calculations. At RHIC energies, 
the non-Gaussian effect was also seen in the experimental HBT correlator, especially in the 
longitudinal direction 2J]. Furthermore, when comparing our fitting result in Fig. [T]with 
that observed in experiment (in Fig. 1 of Ref. [7|), it is seen that the non-Gaussian effect 
is stronger in our calculations than in experiment. This might be due to the omission of 



all potential interactions between partic 
this interpretation was found in Ref. 



es in the current cascade calculations. Support for 
23| , where the consideration of a mean- field potential 
plus Coulomb potential significantly reduced the non-Gaussian effect on correlators of pions 
from HICs at AGS energies. Therefore, both a dynamic treatment of the particle transport 
with a proper equation of state (EoS) for the QGP phase and the hadron phase, and further 
theoretical development of the fitting formalism are equally important for a more precise 
extraction of spatio-temporal information of the source |25[]. 

Besides the non-Gaussian effect, the contribution of the correlation between the emission 



time and position to the HBT radii, especially in the outward direction, has been paid 



more attention in recent years since it closely relates to the stiffness of the EoS of nuclear 



matter especially at the early stage of the whole dynamic process 26N28| . In Fig. |2] (a) 
we show the calculated emission time versus transverse position Xout of pions. The cuts 
\ri\ < 0.8 and 200 < kx < 300 MeV/c are adopted to have the same acceptance as for the 
extraction of the HBT radii. At the same time, since the correlation function and thus 
the HBT radii are mainly sensitive to pairs with small momentum difference, a cut on the 
relative pair momentum \qi\ < 100 MeV/c is applied as well. Further, in order to remove 
the contribution of long-lived resonances, a cut on the freeze-out time (tfreezeout < 199 fm/c) 
is used. It is found that, even in the cascade calculation, there exists a visibly positive 
correlation between the emission time and position. To further analyze the importance of 
the x — t correlation, we quantitatively calculate all three terms in Eq. (|2]) and show them in 
Fig. [2] (b) as a function of kx- As a whole, although the magnitude of the x — t correlation 
term (— 2(/3Ttx) ~ -300 fm^) is as big as that of the emission region term ((x^) ~ 200 fm^), 
the most important contribution to Rq comes from the emission duration term ((/Sl^t^) ~ 
1000 fm^). It implies that both a shorter duration time and a stronger x — t correlation lead 
to a smaller Rq value, which will be further discussed in Fig. [3] specialized for the result of 
HBT radii. Here, it is interesting to see that the direct computation of Rq leads to a value 
of ~ 30 fm which is larger than the value extracted from the Gaussian fit to the correlation 
function by about a factor of three, as was also observed in the AMPT model calculations 
for Au-Au colhsions at RHIC [26 1. 



Fig. [3] shows the kr dependence of the HBT radii Rq, Rs, and Rl, as well as the ratio 
Ro/Rsi extracted from the Gaussian fit to the two-pion correlators. The UrQMD cas- 
cade calculations for central Pb-Pb collisions at LHC energy ^/JNN = 2.76 TeV (lines with 
crosses) and central Au-Au collisions at RHIC energy ^/sNN = 200 GeV (lines with up- 
triangles) are compared to corresponding experimental data by ALICE/LHC (open stars) 
and STAR/RHIC (open diamonds). A strong decrease of the three HBT-radii with k^ is seen 
both in experiments and in the UrQMD calculations for HICs. This implies a substantial 
expansion of the source and is qualitatively captured by the UrQMD dynamics. Following 
experimental results, the calculated HBT radii for Pb-Pb at LHC are found to be larger than 
those for Au-Au at RHIC. The largest increase exists in the longitudinal direction, which is 
also seen by the experiments. Although the comparison of the calculated HBT radii Rl and 
Rs with corresponding data at RHIC is fairly well, it gets worse when going to the higher 
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FIG. 3: (Color online) fcy dependence of pion HBT radii Rq [panel (a)], Rs [(b)], and Rl 
[(c)], as well as the ratio Rq/Rs [(d)], for central {a far < 5%) Pb-Pb collisions at LHC energy 
Y^sjViV = 2.76 TeV. For comparison, parameters for central {cr/aT < 5%) Au-Au collisions at 
RHIC energy ^/sNN = 200 GeV are also shown. Lines with up-triangles and crosses are for model 

calculations while scattered symbols are for experimental data of STAR/RHIC and ALICE/LHC 

Pfl 

collaborations taken from [7|, |2J] . Lines with circles and diamonds show results with an artificially 

decreased emission duration by a factor of tgcaie = 2 and 3, separately, in the analysis of correlation 
function. 

LHC energy. At LHC the calculated Rs values at all kx are found to be slightly smaller 
than data, while Rl and Rq values are larger than data. Together with large calculated Rq, 
the emission time related quantity Rq/Rs is found to be markedly larger than the data. 

From Eq. (|2]) and Fig. [2] (b) it is clear that the value of Rq is strongly dependent on the 
emission duration of the particles. To further investigate the contribution of the emission 
duration to the HBT radii, we artificially decrease it by rescaling the relative time t to the 
"effective source center time" t (=< ti >) by t = tj — t — )■ t' = (tj — t) /tgcaie in the calculation 
of the correlation function at LHC energies. This effectively changes Eq. ([2]) to 



B!^ = {{x - PTt'Y) = {x') - 2 



Xhtx) , (/32t2) 



T^ scale 



+ 



f2 
scale 



(5) 



The results for this calculation are presented as lines with circles {tgcaie = 2) and with 
diamonds {tscaie = 3) in Fig. |3l The artificially decreased emission duration leads to smaller 
Ro values in all kx bins but leaves Rs unchanged, as expected. Overall it results in an 
improved agreement with the data of Rq/Rs ratio. From Fig. 12] it is also found that R^ is 
overestimated at LHC. Since Rl is mainly related to the lifetime of the source, it implies 
that this lifetime is also overestimated by UrQMD. Other calculations in |7|, l29| show that 
UrQMD overestimates the source lifetime by a factor of ~ 2 — 3 when compared to LHC 
data. The overestimation of both Rq and Rl can be attributed to the known fact that 
the pressure in the early stage is not strong enough in the cascade model calculations. A 
higher pressure would lead to a more explosive expansion, a stronger phase-space correlation, 
and a faster decoupling of the system, thus leading to smaller regions of homogeneity. For 



hydrodynamics hybrid approach |14j . Il8l . 



more discussion we refer the reader t o | 25l. 27 1. With the improved integrated Boltzmann + 



31| , where various EoS of nuclear matter during 



the hydrodynamic evolution may be treated consistently and a decoupling supplemented by 
realistic 3d hypersurfaces we hope to get a satisfactory solution in the near future. 

IV. SUMMARY 

To summarize, the two-pion HBT correlations (in the LCMS system) for central Pb-Pb 
collisions at the LHC energy ^/snn = 2.76 TeV are calculated for the first time with the 
microscopic transport model UrQMD. The non-Gaussian effect is seen especially in both 
longitudinal and outward directions. Both the transverse momentum kj- dependence and 
the beam energy (from RHIC to LHC) dependence of the HBT radii Rq, Rs, and Rl, 
extracted from a three dimensional Gaussian fit to the correlator, exhibit qualitatively the 
same behaviour as found in the experiments. However, the calculated Rq/Rs ratios at all 
kr bins are found to be larger than in the data, both at RHIC & LHC. We traced this 
finding back to the explosive dynamics of the fireball at LHC which results in both a shorter 
emission duration and a stronger time-space correlation than modeled here. 
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